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New publications (http://shop.bk-teachware.com) 

[1] H-J Kayser SR-52: Lehr- und Lernvideos zur Mathematik mit Derive, bk-teachware 

[2] H-J Kayser. SR-53: Projekte im Mathematikunterricht mit Derive, bk-teachware 

[3] Bernhard Kutzler. SR-B: Die Sprache der Zahlen, bk-teachware 


Announcement of a new publication: 

Gutierrez, A., Boero, P. (Eds.). (2006). Handbook of Research on the Psychology of Mathematics 
Education. Rotterdam, The Netherlands: Sense Publishers. 

This volume is a compilation of the research produced by the International Group for the Psychology 
of Mathematics Education (PME) since its creation, 30 years ago. It has been written to become an 
essential reference for Mathematics Education research in coming years. 

The chapters offer summaries and synthesis of the research produced by the PME Group, pre- 
sented to let the readers grasp the evolution of paradigms, questions, methodologies and most relevant 
research results during last 30 years. They also include extensive lists of references. The chapters raise 
also the main current research questions and suggest directions for future research. 

The volume is the result of the effort of 30 authors and 26 reviewers. Most of them are recog- 
nized leading international researchers, members of the PME Group, with great expertise on the topic 
of their chapter. It includes 15 chapters, divided into five sections, devoted to the main research do- 
mains of interest to the PME Group. The first to third sections summarize cognitively oriented re- 
search on learning and teaching specific content areas (algebra, geometry and measurement, numerical 
thinking), transversal areas (ATM, proof, visualization, young children’s mathematical thinking), and 
based on technology rich environments (use of technology for teaching and learning algebra, calculus, 
and geometry). The fourth section is devoted to the research on social, affective, cultural and cognitive 
aspects of Mathematics Education (affectivity, constructivism, equity, socio-cultural practices). Fi- 
nally, the fifth section includes two chapters summarizing the PME research on teachers training and 
professional life of mathematics teachers. 

This handbook shall be of interest to both experienced researchers and doctoral students needing 
detailed synthesis of the advances and future directions of research in Mathematics Education, and 
also to mathematics teacher trainers who need to have a comprehensive reference as background for 
their courses on Mathematics Education. You can get more information of the book, and download a 
nonprintable PDF version of it, from Sense Publishers’ web page at 
http://www.sensepublishers.com/books/otherbooks/90-77874-19-4.htm . 


Josef Bohm, here is a link to tables of elliptic integrals generated with DERIVE. 

http://www.getnet.net/-cherrv/derive/index.htnnl 

Regards, Jim FitzSimons 


Download all DNL-DERIVE- and Tl-files from 

http : / / www . austromath . ac . at/ dug/ 
http : / / www . derive-europe . com/ support . asp?dug 
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Dear DUG Members, 

I am happy that I can publish DNL#62 some days 
before DES-2006 in Dresden will start. This is a 
very special Newsletter and it deserves its name 
“Bulletin of the DERIVE User Group”, because all 
contributions developed from User requests and the 
respective answers. 

The extended exchange of emails concerning 
COS(7i/17) induced “Mr Titbits” Johann 
Wiesenbauer to bring Derive on its calculation and 
manipulation boundaries in his Titbits #32. Jan 
Vermeylen sent a challenge for applying a slider 
bar to animate the motion of a particle in a mag- 
netic field given by some thousand positions in 
space. An international cooperation between Swe- 
den, England and Austria found a satisfying solu- 
tion. Many thanks to David Sjostrand and Peter 
Schofield. Finally, Heinz Rainer Geyer reminded 
me on my career as a teacher with his problem 
finding a generalised rule to present the frequency 
distribution of the sum of n rolled dice. Here again 
Johann Wiesenbauer found the shortest and most 
clever solution. 



Unfortunately we are very poor in Tl-contributions. 
Wolfgang Propper announced a Tl-paper based on 
a lecture which he gave in Vienna some months 
ago. Many thanks in advance. Ed like to demon- 
strate that in many cases it would not be too diffi- 
cult to transfer the problems and their solutions 
onto the handheld (dice problem, random walk). 


There are other interesting problems (together with 
their solutions) remaining on my desk: how to load 
an external data set to DERIVE and assigning it 
simultaneously to a variable name. Klaus Rohe 
raised this question and he promised a contribution 
on the “2 and 3 Body Problem” applying the pro- 
posed solution. Another request was on applying 
slider bars to Peter Schofield sophisticated 
“ARROW”-tool (a user contributed DERIVE utility 
file) and there was a also question from a Univer- 
sity lecturer on a special representation of Taylor 
series. 

By the way, do you now how to “paint” in and 
from Derive? A student of Tania Roller found this 
out by experimenting with the program. 

The two pictures are from a workshop with 50 (!!!) 
students on the Vienna International School. I was 
invited by Marlene Torres-Skoumal to hold a 
“Modelling Session” with the 10-graders of this 
highly reputated institution. I felt very fine among 
students from so many countries and we agreed to 
have another workshop in the next term. 



Finally Ed like to wish you all a wonderful summer 
and I hope to meet many of you in Dresden. 



Invitation 

DERIVE User Group Meeting 

DES-TIME 2006 in Dresden 
Saturday 22 July 2006, afternoon 

Please follow the announcement given on the Conference 


EDITORIAL 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS- 77 User Group. 
It is published at least four times a year 
with contents of 44 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE 
and the 77-89 /92/ Titanium! Voyage 200 as 
well as to create a group to discuss the 
possibilities of new methodical and didac- 
tical manners in teaching mathematics. 

As many of the DERIVE Users are also 
using the CAS-TIs the DNL tries to com- 
bine the applications of these modem tech- 
nologies. 


Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone/FAX: 43-(0)2275/8207 
e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS-77 Newsletter will be. 


Next issue: September 2006 

Deadline 3 1 August 2006 


Preview: Contributions waiting to be published 

Two Stage Least Squares, M. R. Phillips, USA 
Some simulations of Random Experiments, J. Bohm, AUT 
Wonderful World of Pedal Curves, J. Bohm 
Another Task for End Examination, J. Lechner, AUT 
Tools for 3D-Problems, P. Luke-Rosendahl, GER 
ANOVA with DERIVE & Tl, M. R. Phillips, USA 
Financial Mathematics 4, M. R. Phillips 
Hill-Encription, J. Bohm 

Farey Sequences on the Tl, M. Lesmes-Acosta, COL 
Simulating a Graphing Calculator in DERIVE, J. Bohm 
Henon & Co, J. Bohm 

Challenges from Fermat, Bj. Felsager, DEN 
Are all Bodies falling equally fast, J. Lechner 
Modelling Traffic Density, Th. Himmelbauer, AUT 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 
An Interesting Problem with a Triangle, Steiner Point, P. Luke-Rosendahl, GER 
Mathematics and Design, Hubert Weller, GER 
Diophantine Polynomials, Duncan E. McDougall, Canada 
Contour Plots and Implicit Plots, Peter Schofield, UK 

and Setif, FRA; Vermeylen, BEL; Leinbach, USA; Koller, AUT; Baumann, GER; 
Keunecke, GER and others 

Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 
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From cos( 7 i/ 17 ) to acot(x) 

An interesting exchange of emails 


Starting the discussion: 

How does DERIVE find: 


14 May 2006 


COS 


JUOS-J17 + 170) + 3-47 + 17) 7(34 - 2-../17) 

+ 


17 


S 


16 


Regards, Jim FitzSimons 


-J17 1 

+ 

16 16 


Hello Jim, 


15 May 2006 


A quick way to get an approximative value of cos(pi/17) is to first simplify 

APPROX (NSOLVE (CHEBYCHEV_T (17, x) = -1, x, 0.9, 1), 30) 

and then to approximate the output, which should lead to the approximation 

X = 0.9829730996 for COS(ti/17). 

In order to get your radical expression I'm afraid you can't avoid following in Gauss tracks, 
when he computed it as a 19 year old. 

Unfortunately it is too long to describe it a few sentences, but I have given a detailed descrip- 
tion of his approach in my Titbits series in the Derive Newletter#20. 

Hope this helps. 

Cheers, 

Johann (Wiesenbauer) 

Jim asked for a copy of Johann 's Titbits from DNL#20 and I sent him this contribution. Jim joined the 
DERIVE User Group and we give him a warm welcome in our group, Josef 


Hello Jim, 


15 May 2006 


solve the equation z A 17 = —1 in Polarform! 

You get z = cos(pi/17 + 2/17*k*PI) + i sin(PI/17 + 2/17*k*PI). For k=0 the Real-part of the 
solution leads to cos(PI/17). 

I don't know if Derive does it this way! But I would! 


Greetings, 

Manuel 
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Hello Manuel et al., 


21 May 2006 


Sorry for the delayed answer, but it took me some time to "fill in the details" as to your pro- 
posal of a computation of cos(7t/17) using Derive. It is essentially following Gauss' solution 
and the outline I gave of it in my Titbits(6) in the DNL #20. (For Jim: This is the Newletter of 
the Derive User Group, some issues of which can be found online using the keywords dug 
derive in a Google search.) 

If you (or anybody else reading this) can find a shortcut somewhere, I would be very much 
interested in it. 


Cheers, 

Johann 

File cospil7 . dfw 


Computation of cos(pi/17) using Derive 


((c) Johann Wiesenbauer, Vienna University of Technology, 2006) 


We start with the observation that cos(pi/17) is a solution of the polynomial equation 
chebychev_T(17, x)+l =0 (*) 

due to the defining relation chebychev_T(n,cos <|>) = cos(n<|>). In the following we use 
Dickson polynomials D(n,x) though, which are closely related to Chebychev polynomials 
of the first kind. More precisely, the relation 

Dickson(n,x) = 2*chebychev_T(n, x/2) 

holds and therefore 2cos(pi/17) is a solution of the polynomial equation 
Dickson(17,x)+2=0 (**) 

The polynomials Dickson(n,x) can be also computed in Derive using the definition 

#1: Dickson (n, x) := V_M0D(n, x, 1) 

Hence the LHS of (**) is given by the polynomial 

#2: pi := Dickson(17, x) + 2 

17 15 13 11 9 7 

#3: pi := x - 17 -x + 119-x - 442-x + 935-x - 1122-x 

5 3 

+ 714 -x - 204-x + 17 -x + 2 


D 
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Factoring over the field Q of rational numbers yields 

#4: FACTOR (pi , Rational, x) 


8765432 2 

#5: (x + 2) • (x - x - 7-x + 6-x + 15-x - 10-x - 10-x + 4-x + 1) 


In fact, this is the factorization of pi into irreducible factors and the 8-degree factor 
p2 is the minimal polynomial of 2cos(pi/17) over Q (as well as of its conjugates 
2cos(3pi/17),2cos(5pi/17) 2cos(15pi/17)). 

8 7 6 5 4 3 2 

#6: p2:=x -x - 7*x + 6-x +15*x - 10*x - 10*x + 4*x +1 


In order to show this we define the vector cosines consisting of all those cosine values. 
For example, 2cos(pi/17) = z - z A 16 with z = exp(#i pi/17). (Note that we never define z 
explicitly, but only use (z~17+l)/(z+l)=0 in the following!) 


#7: 


reduce(u) := REMAINDER 


u, 


17 

z + 1 


z + 1 J 


k 17 - k 

#8: cosines := VECT0R(z - z , k, 1, 16, 2) 

T 16 3 14 5 12 7 10 9 

#9: cosines :=Lz-z ,z - z ,z - z ,z - z ,z - 

8 11 6 13 4 15 21 

z,z -z,z -z,z -zj 

VECT0R( reduce (SUBST(p2 , x, cosines )), k, 1, 8) = [0, 0, 
#10: k 


0 , 0 , 0 , 0 , 0 , 0 ] 


By the way, Derive thinks that p2 cannot be factored using radicals, i.e. 

8 7 6 5 4 

#11: FACT0R(p2, Radical, x)=x -x - 7*x +6*x +15*x - 

3 2 

10*x - 10*x + 4*x + 1 

In the following we'll show that this is wrong. In fact, it will turn out that p2 can be 
split up into two polynomials of degree 4 in the ring of integers of Q[sqrt(17)]. 

In order to show this we first observe that the sum of all cosine values must be the 
negative coefficient of x~7 of p2, i.e. it must be 1. 
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# 12 : 


reduced(cosines)) = 1 


Now we face the following problem: Is it possible to break up this sum into two partial 
sums with 4 summands each such that their product is in Q? Do you think that Derive 
can help us out here? 

SELECT(RATIONAL?(reduce(I(cosines )-X(cosines 

#13: SORT(s) SORT(~ 

))), s, P0WER_SET ({1 , 8}, 4)) 

{1 8} \ s) 

#14: {{1, 5, 7, 8}, {2, 3, 4, 6}} 


Wow! There is exactly one such partition and Derive found it in splits of a second! In 
fact, we have 


[ si := Hcosines ), s2 := Kcosines )] 

[1, 5, 7, 8] [2, 3, 4, 6] J 

#16: reduce(sl*s2) = -4 


which means that the real values corresponding to si and s2 must be the among the 
roots of 


#17: S0LUTI0NS([u + v = 1, u-v = -4], [u, v]) 


#18: 


Jl 7 1 1 Jl7 

2 2 2 2 

1 417 4 17 1 

— - + — 

2 2 2 2 


There is a small problem though: We don't know the correspondence between solutions 
u,v and sl,s2. One way to find this out is to use their numerical values: 


#19: APPROX 


417 1 

+ 

2 2 

1 417 

2 2 


1 417 

2 2 

417 1 

— + — 

2 2 


2.561552812 -1.561552812 ' 

-1.561552812 2.561552812 . 


#20: APPROX 

2 

r 

cos 

' C2-k - 1)-7T " 

N 

, k, {1, 5, 7, 8} 




c 

, 17 

y 


r 

' C2-k - 1)-7T ' 

N - 

cos 


, k, {2, 3, 4, 6} 


, 17 

y . 


# 21 : 


[-1.561552812, 2.561552812] 
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Hence, as from these numerical values, we have got the following equations 
l 


# 22 : 


si = 


yi 7 417 i 

, s 2 — + — 

2 2 2 J 


Furthermore, the partition {{1, 5, 7, 8}, {2, 3, 4, 6}} of (1,2, 3, 4, 5, 6 ,7, 8} tells us exactly 
the grouping of linear factors of p2 for the factorization of p2 mentioned above. In the 
first place, it's easy to get the approximative representation of the factors by 


#23: APPROX 

s 

EXPAND 

n 

x - 2- COS 

' (2 • k - 1).7T ' 

, k, {1, 5, 7, 8} 






, 17 

) 

// 


4 3 2 

#24: x + 1. 561552812- x - 3 . 561552812 -x - 6.123105625-x - 1 


#25: APPROX 

r 

EXPAND 

r 

n 

r 

x - 2- COS 

' (2-k - 1)-7T " 

, k, {2, 3, 4, 6} 






, 17 

) 



4 3 2 

#26: x - 2. 561552812- x + 0.5615528128-x + 2.123105625-x - 1 


But what about their exact representations? Here we make use of the fact that the 
exact coefficients are linear combinations of si and s2. (For the cognoscenti, {sl,s2} is a 
basis for the field extension Q(sqrt(17)) over Q.) Hence, we could make a list of "small" 
linear combinations and simply compare its entries in the first colums with numerical 
values of the coefficients of our polynomials above, e.g. 


#27: 



417 417 1 ' 

1 + 

2 2 2 . 


#28: 


APPR0X(S0RT (VECT0R([ s • x , x], x, 
' -12.36931687 
-10.80776406 
-9.807764064 
-9.246211251 
-8.246211251 
-7.684658438 
-7.246211251 
-6.684658438 
-6.123105625 


{-3 3} ))) 

[3, -3] ' 

[2, -3] 

[3, -2] 

[1, -3] 

[ 2 , - 2 ] 

[0, -3] 

[3, -1] 

[ 1 , - 2 ] 

[-1, -3] 


[ 2 , - 1 ] 


-5.684658438 


- 5.123105625 

[ 0 , - 2 ] 

- 4.684658438 

[ 3 , 0 ] 

- 4.561552812 

l 1 

1 

hJ 

1 

UJ 

i i 

- 4.123105625 

[ 1 , - 1 ] 

- 3.561552812 

[- 1 , - 2 ] 

- 3.123105625 

[ 2 , 0 ] 

-3 

i 1 

1 

UJ 

1 

UJ 
1 1 

- 2.561552812 

1 — 1 

0 

1 

M 

i i 

- 2.123105625 

[ 3 , 1 ] 

-2 

i 1 

1 

l\J 

1 

l\J 

i i 

- 1.561552812 

[ 1 , 0 ] 

-1 

[- 1 , - 1 ] 

- 0.5615528128 

[ 2 , 1 ] 

- 0.4384471871 

l 1 

1 

UJ 

1 

l\J 

i i 

0 

[ 0 , 0 ] 

0.4384471871 

[ 3 , 2 ] 

0.5615528128 

[- 2 , - 1 ] 

1 

[ 1 , 1 ] 

1.561552812 

[- 1 , 0 ] 

2 

[ 2 , 2 ] 

2.123105625 

i 1 

1 

UJ 

1 

I— 1 

i i 

2.561552812 

[ 0 , 1 ] 

3 

[ 3 , 3 ] 

3.123105625 

[- 2 , 0 ] 

3.561552812 

[ 1 , 2 ] 

4.123105625 

[- 1 , 1 ] 

4.561552812 

[ 2 , 3 ] 

4.684658438 

[- 3 , 0 ] 

5.123105625 

[ 0 , 2 ] 

5.684658438 

[- 2 , 1 ] 

6.123105625 

[ 1 , 3 ] 

6.684658438 

i 1 

1 

M 

hJ 

i i 
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7.246211251 

[-3, 

1] 

7.684658438 

[ 0 , 

3] 

8.246211251 

[-2, 

2] 

9.246211251 

[-1, 

3] 

9.807764064 

[-3, 

2] 

10.80776406 

[-2, 

3] 

. 12.36931687 

[-3, 

3] . 


Looking up this table we get the following "exact" factors of p2: 

4 3 2 

#30: p21 := x + ( s - [ —1 , 0])-x + (s*[-l, -2])-x + (s*[-l, 


-3])*x + s • [ -1 , -1] 


#31: 


#32: 


4 3 

' 717 

1 ' 

2 

' 7l 7 3 > 

x + x • 

, 2 

2 , 

- x • 

,2 2 > 


- x • (Jl7 + 2) - 1 


4 3 2 

p22 := x + (s * [ 0 , -l])-x + (s • [-2 , -l])-x + (s - [-3 , 


-l])-x + s • [ -1 , -1] 


4 3 

#33: p22 := x - x • 


'717 1 ' 

2 

717 

3 ' 

+ 

2 2 J 

+ x • 

, 2 

2 , 


+ x-G/17 - 2) - 1 


Indeed, Derive can verify this factorization now (though only after one of those stran- 
ge "double simplifications"!) 

#34: p2 - p21*p22 


#35: 

#36: 

Let's continue by setting 

#37: p3 := p21 


4 3 

#38: p3 := x + x • 


x .((^17 - 2). (717 + 2) - 13) 
0 


717 

1 ' 

2 

'717 3 

, 2 

2 , 

< 

i 

,2 2 > 


x • (7l7 + 2) - 1 


Now by the same reasoning as above, we know that the sum of its zeros, namely 

#39: Kcosines 


16 15 13 9 8 4 2 

) = — z +z +z +z - z - z - z + z 


[1, 5, 7, 8] 
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is an element of Q(sqrt(17)), but how can we break up this sum into two sums of equal 
size such that their product is in Q(sqrt(17)) too? The following computation shows that 
the partition {{1,7}, {5, 8}} works: 


reduce(I(cosines )-X(cosines )) = -1 

#40: [1, 7] [5, 8] 


It should be clear by now how the story goes on. 


/ 

1 

717 

N 

#41: SOLUTIONS 

U + V = - 

, U • V = -1 

, [u, v] 

s 

2 

2 

y 


#42: 


7(34 - 2-717) 717 1 

- + 

4 4 4 


7(34 - 2-717) 


7(34 - 2-717) 717 1 7(34 - 2-717) 

4 4 4 4 


717 1 

+ 

4 4 

Vi 7 l 

4 4 


#43 : APPROX 


V(34 - 2-Vl 7) Vl7 1 

4 4 4 

V(34 - 2-Vl 7) Vl7 1 

- + 

4 4 4 


V(34 - 2 • Vl7) Vl 7 1 

— + — 

4 4 4 

7(34 - 2-717) 717 1 

4 4 4 


0.4879283649 -2.049481177 

-2.049481177 0.4879283649 


#44: APPROX 



r r 

2-1 

cos 

. 

( V 


(2 • k - 1)-7T "1 


17 


, k, {1, 7} 


2-1 


COS 


(2 • k - 1) -7T 


17 


, k, {5, 8} 


= [0.4879283649, -2.049481177] 
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From the approximations above we can conclude that 


#45 


[tl := Hcosines ), t2 := Hcosines )] 

: L [1, 7] [5, 8] J 


#46: 


7(34 - 2-717) 717 1 

tl = - + , t2 = - 

4 4 4 


7(34 - 2-717) 717 1 

- + 

4 4 4 


Finally we have 

14 12 5 3 

#47: reduce(cosi nes -cosines )=z +z -z-z 

1 7 

11 10 7 6 

#48: reduce(cosines -cosines )=-z +z -z +z 

5 8 

which can also be written as 



reduce(cosi nes -cosines ) = 

- Hcosines 

) 

#49: 

1 7 

[2, 

3] 


reduce(cosi nes -cosines ) = 

- Hcosines 

) 

#50: 

5 8 

[4, 

6] 


The sums on the RHS of these equation can be computed as follows 


frl := Hcosines ), r2 := Hcosines )] 

#51: L [2, 3] [4, 6] J 

717 1 

#52: rl + r2 = + 

2 2 


#53: reduce(rl- r2) = -1 


The rest is supposed to need no further comment. 


T 717 i 


U + V = + , U • V = -1 

l 1 

C 

< 
i i 

U 2 2 

y 


#54: SOLUTIONS 
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#61: 


#62: 


1.965946199 -1.478017834 

-1.478017834 1.965946199 J 


APPROX 


r 

r 71 V 

2- COS 



, 17 1 


= 1.965946199 


#63 : COS 


7(7(38-717 + 170) + 3-717 + 17) 


17 J 

7(34 - 2-717) 7l 7 


8 


16 16 16 

And what does Derive think of this final equation? Let's conclude with 

7(7(38-717 + 170) + 3-717 + 17) 


#64: APPROX 


COS 


L 17 

7(34 - 2-717) 717 


16 


16 


16 


8 


= true 


21 May 2006 


Johann, thank you for your help. I have a copy of TITBITS6. 

Here are links to my first try. 

http : / / www2 . pvc . maricopa . edu/~f it z simons /cosine . doc 
http : / / www2 . pvc . maricopa . edu/~f it z simons/ cosine . mth 


The college web site is down now, but it will be back up soon. 
Regards, Jim 


Hello all, 


22 May 2006 


Is there any prime p>17 such that COS(7i/p) or SIN(7i/p) is expressible in terms of radicals? 
Aloha, 

Albert D. Rich 


Well, there are two more Fermat primes known: 257 and 65537. 
Cheers, 

David W. Cantrell 


22 May 2006 
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22 May 2006 

Hello Albert, 

As David has already mentioned such a prime must be of Fermat type, hence p=257 and 
p=65537 are certainly further solutions. It is not even known whether are finitely many Fer- 
mat primes, although this is probably true for heuristic reasons. 

This also gives me the opprortunity to correct a mistake in my presentation (near line #40) in 
the attached file below, which went undetected due to a strange coincidence and didn't affect 
the computations thereafter. 

Cheers, Johann 

( Comment of the publisher: the printed file cospi17.dfw is the latest version.) 

22 May 2006 

Hello all, 

Thank you for your positive response to my question. 

Why does a Fermat prime p of the form2 A 2 A n+l, where n is a nonnegative integer, guarantee that 
COS(pi/p) is expressible in terms of radicals? 

Has anyone ever bothered to determine the radical representation of COS(7i/257)? It must be huge... 
Aloha, Albert 


23 May 2006 

Hello Albert, 

First a word of caution: "Expressible in terms of radicals" means here that a term can be 
achieved starting with rational numbers and applying the operations +,-,*,/, SQRT() finitely 
many times in any order. In particular, radicals other than square roots are not allowed. Hav- 
ing said this what follows is an outline of the proof: 

Suppose that u = cos(pi/p) is expressible in terms of radicals for some odd prime p. Then the 
minimal polynomial m(x) of cos(pi/p) over the field Q of rational numbers is given by 

m(x) = (x A p+1)/(x+1) 

as can be seen by applying Eisenstein's criterion to the polynomial m(x-1 ). Hence, the alge- 
braic field extension Q(u) is of degree p-1 . On the other hand, each operation when "build- 
ing" the term for u corresponds to a field extension of degree 1 or 2 and multiplying all those 
degrees together we get a power 2 A k, which must be equal to p-1 . Hence, p is a prime of the 
from 2 A k+1 , which implies that k = 2 A n for some n and p is a Fermat prime. Conversely, if p is 
a Fermat prime F_n, then the Galois group of the field extension Q(u):Q is cyclic and of order 
2 A 2 A n. This means that this extension can be broken up into a series of extensions of degree 
2, leading to a radical expression for u. 

Hope this helps. 

Cheers, Johann 

PS. As for the radical expression of cos(pi/257), it is huge, indeed. I'll try to achieve it in De- 
rive, but this may take some time... 
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23 May 2006 

Thank you, Johann, for making that explicit. However, that raises a question I had, but did not men- 
tion, when I first replied to Albert: 


He had asked if there were any primes p > 17 for which cos(pi/p) is expressible in terms of radicals. I 
knew there were, and replied, merely based on my knowledge of constructibility of regular polygons, 
that there are two Fermat primes >17 known. I realized that the radicals required in those two cases 
would be just square roots. 


You said "As David has already mentioned such a prime must be of Fermat type”, but I didn't say 
quite that. In fact, I was careful to avoid saying that such a prime must be of Fermat type because I 
thought that Albert was asking about expressibility in terms of general radicals, not necessarily square 
roots. 


What I don't know is whether allowing other radicals (cube roots, fifth roots, etc., and thus taking us 
out of the realm of constructible regular polygons) would allow us to express cos(pi/p) in terms of 
radicals for other primes p. Do you know the answer to that? 

BTW, Albert later asked if anyone had ever determined the radical expression for COS(tc/ 257). Those 
interested in that question might like to see 


< http://mathworld.wolfram.com/257-gon.html > and 
< http://mathworld.wolfram.com/65537-gon.html >. 

Regards, 

David 


23 May 2006 

Ah, silly me! I just hadn't thought about the question much beforehand, but I now see that the answer 
is YES. For example, COS(7i/7) can be expressed in terms of radicals (specifically, square and cube 
roots). 

David 


23 May 2006 

For anyone who might be interested, I've now gotten COS(7i/7) in a nice form. Letting 
u = (7/2 (-1 + 3 sqrt(-3))) A (l/3), we have COS(tt/ 7) = 1/6 (u + 1 + 7/u). 

David 


23 May 2006 

David, I know you can express COS(7i/7) in cube roots. COS(tt/ 7) can not be expressed in square 
roots. 

Regards, Jim FitzSimons 


23 May 2006 

Thanks Johann and David for your insightful responses to my questions. 

If u := (7/2 (-1 + 3 V(-3))) A (l/3), Derive simplifies 


COS(tt/ 7) = 1/6 (u + 1 + 7/u) 
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to 

COS(tt/7) = V(7)-COS(ACOT(-V(3)/9)/3)/3 + 1/6 

This leads to the potential simplification rule 

COS(ACOT(-V(3)/9)/3) -► V(7)-(6-COS(tc/7) - 1)/14 

That leads me to wonder if there is a generalization of this rule that would simplify expressions of the 
form COS(ACOT(n)/3) and COS(ATAN(n)/3) where n is an algebraic number to expressions free of 
inverse trig functions. Any thoughts on this? 

Aloha, 

Albert 


23 May 2006 

It all depends, I think, on whether we require that the radicals be real. If not, then we can say, for ex- 
ample, 

CAS(ACOT(x)/r) = 1/2 (u + 1/u) 

where u = ((x + #i)/ V(x A 2 + l)) A (l/r). 

Is that the sort of thing you wanted? 

Cheers, 

David 


24 May 2006 

Thanks for pointing out the nice identity for COS(ACOT(x)/r). I wasn't aware of it. 

Unfortunately, the answer to your question is 'no'. For it to be a transformation rule Derive can use, 
the expression equivalent to COS(ACOT(n)/3) must simpler in the sense that it consists only of real 
radicals and COS(7i/m) where m is an integer. I doubt there exists such an expression for arbitrary 
algebraic numbers n; however as the example below shows, there are at least some n for which there is 
a simpler equivalent. 

So my question is: What is the class of algebraic numbers for which there is a simpler equivalent 
AND what is that equivalent? 

Aloha, 

Albert 


24 May 2006 

Hallo Albert, 

There must be something wrong with this identity because a simple 
numerical observation leads to the following contradiction within 
DERIVE : 

(1) COS (ACOT (-SQRT (3) /9) /3 is approximately 0.8326204335 

(2) COS (ACOT (SQRT (3) /9) /3 is approximately 0.8959532196 
These are the results due to Derive 5.06 and 6.01. 
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But both results must be equal, because arccot (-x) =-arccot (x) and 
cos (-x)=cos (x) . So there appears a numerical discrepancy. 

On the other hand SQRT ( 7 ) * ( 6*C0S (pi/7 ) - 1)/14 evaluates numeri- 

cally to 0.8326204335. 

Other CAS and numerical calculators suggest that both, 

COS (ACOT (-SQRT (3) /9) /3) and COS (ACOT ( SQRT ( 3 ) / 9 ) /3 ) , should evaluate 
numerically to 0.8959532196, and that SQRT ( 7 ) * ( 6*COS (pi/7 ) - 1)/14 
should evaluate numerically to 0.8326204335. 


So I conclude at a first glance that the identity above may be 
false . 


"Aloha", Stefan Welke 


24 May 2006 

Albert, I searched for solutions for COS(x/3) and I have not found any yet. 

Regards, Jim FitzSimons 


24 May 2006 

Hi Jim, 

I am concerned about the special case when your x is of the form ACOT(n) where n is an algebraic 
number. For example, as shown below COS(ACOT(n)/3) does simplify for n = -SQRT(3)/9. 

The question is: What is the class of n for which COS(ACOT(n)/3) simplifies AND what exactly is 
the simplification rule? 

Aloha, Albert 


24 May 2006 

Hello Stefan, 

Your assertion that ARCCOT(-x) = - ARCCOT(x) is incorrect. As plotting quickly shows, the cor- 
rect identity is 

ARCCOT(-x) = ti- ARCCOT( x) 

Therefore, Derive's approximations of the below expressions are both correct, and the identity I gave 
is also correct. 


Aloha, Albert 


24 May 2006 

Albert, 

Let me try to nip a potential argument in the bud. 

You told Stefan that his ’’assertion that ARCCOT(-x) = -ARCCOT(x) is incorrect.” You are undenia- 
bly correct, but only because this is DERIVE-NEWS. 

Let us not forget that there are two common ways of defining ARCCOT(x), each with its own advan- 
tages and disadvantages. Derive uses one of those ways. 

Done the other way, ARCCOT(-x) = -ARCCOT(x) holds for all nonzero real x. 


David 
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InputMode:=Word 


25 May 2006 


"The question is: What is the class of n for which C0S(AC0T(n)/3) simplifies 

AND what exactly is the simplification rule?" 


m=[l, 2, 3, 4, 5, 6, 8, 10, 15, 16, 17, 20,". . ."] 


k: epsi lonlnteger 
C0S(AC0T(n)/3) 

AC0T(n)/3=k*pi /m 
;Solve(#6,n) 

SOLVE (ACOT(n) /3=k*pi /m , n) 

;Simp(Solve(#6,n)) 

n=IF((3*k-m)/m<0 AND (6*k-m)/m>-l,C0T(3*pi*k/m)) 


Regards, Jim FitzSimons 


25 May 2006 

n - -SQRT(3)/9. 

Gives COS(7i/7) which can not be expressed as with only sqrt radicals. 

Regards, Jim FitzSimons 


25 May 2006 

Yes, I should have been more precise by adding the phrase "in Derive, at least" to the end of the sen- 
tence you quote below. Thanks for nipping this argument in the bud! 

However, now Fm curious as to what percentage of DeriveNews readers prefer the Derive definition 
of ARCCOT(x), and what percentage don't. For simplicity sake, let's stick to real-valued x. 

Derive defines ARCCOT(x) as tt/2-ARCTAN(x). Therefore ARCCOT(-x) equals tt-ARCCOT(x). 
How is the alternative definition of ARCCOT(x) defined in terms of ARCTAN(x) so that 
ARCCOT(-x) equals -ARCCOT(x)? 

Aloha, 

Albert 


Hello Albert, 


25 May 2006 


sorry, I should had done a plot of ACOT in DERIVE first. What I had 
to learn as a mathematician and as a teacher is that most CAS and 
calculators define arccot (x) as arctan (1/x) which is technically 
very simple because you only have to implement arctan (or arccot) . 
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One reason may be the following argument: y=arccot (x) is equivalent 
to x=cot(y). We know that tan (y) =l/cot (y) =l/x, so 
y=arctan ( tan (y) ) =arctan ( 1/x) . We finally arrive at arc- 
cot (y) =arctan ( 1/x) . 


In pre-computer times I learned at school and university that the 
real cotangent function must be first restricted to the open inter- 
val ]0,pi[ where it is strictly monotone. Now the restricted cotan- 
gent function has an inverse arccot, which is defined on ] -inf, inf [ 
with arccot (] -inf , inf [ ) =]0,pi[. This explains why arc- 
cot (y) =arctan ( 1/x) is mathematically incorrect. The simple reasoning 
above ignores the fact that the inverse functions arctan and arccot 
are defined only for restrictions of the the tan and cot functions. 
By the way, it is very difficult to teach to average students that 
they have to interpret the results of their calculators with these 
facts in mind. 


To answer your question, Albert: I am very comfortable with the way, 
DERIVE defines the arccot function because it is mathematically cor- 
rect for real numbers . 

What I learned once again is the importance to know, how the par- 
ticular program you are working with, defines its functions. This 
knowledge is indispensable to obtain correct and reliable results. 

Aloha, 

Stefan Welke 


25 May 2006 

Albert, ARCCOT(x) as nil - ARCTAN(x) is the correct definition since it gives ARCCOT(x) the 
correct range. 

I teach trigonometry and in my books the domain of arccot(x) is all real and the range is 0 to n. 
Regards, Jim FitzSimons 


Hello Albert, 


25 May 2006 


M For simplicity sake” indeed! Perhaps you don't realize it, but choosing to "stick to real-valued x" is 
highly prejudicial in this discussion. AFAIK, any sane person would choose Derive's definition if we 
are to deal only with real values. And that is surely why Jim thinks that Derive's is "the correct defini- 
tion". But, for better or worse, there simply is no single "correct" definition if we go beyond the reals. 
Out of the infinitude of correct ways that we could restrict the multivalued inverse cotangent relation 
in order to get a single valued function, there are two that stand out as most convenient. Alas, we 
mathematicians have no international standardizing body (or at least none that I'm aware of!), and so 
two common conventions persist. 

I should also note perhaps, for anyone not familiar with it, that the only trig functions which have just 
a single convention for their inverse function are sine, cosine and tangent. The other three all have two 
distinct conventions. For this reason, I tend to avoid using arccot, arcsec and arccsc. This is always 
easily enough done. After all, anything that can be expressed using arccot can also be expressed using 
arctan instead. And done the latter way, there can be no possible ambiguity in the expression. 
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Hello Stefan, 

I agree with almost everything you said. 

> sorry, I should had done a plot of ACOT in DERIVE first. What I had to learn as a mathematician 

> and as a teacher is that most CAS and calculators define arccot(x) as arctan(l/x) which is techni 

> cally very simple because you only have to implement arctan (or arccot). 

First, I'm not sure that most define it that way. I haven't done a poll. But I suspect that you're correct. 

Second, in order to have a complete definition, one cannot simply say "define arccot(x) as arctan(l/x)" 
because that leaves arccot(O) undefined. It gives us no idea whether we should choose arccot(O) = 7i/2 
or arccot(O) = -n/2. 

BTW, Stefan (assuming, perhaps incorrectly, that you're German), does DIN specify a definition of 
arccot as being the only "correct" one? 

Regards, 

David 

25 May 206 

On Thu, 25 May 2006 15:31:15 +0100, David W. Cantrell wrote: 

>Alas, we mathematicians have no international standardizing body (or at 
>least none that I'm aware of!), and so two common conventions persist. 

ISO/IEC 10967-2:2001 (Language Independent Arithmetic, part 2: 

Elementary numerical functions; (LIA-2)), section 5.3.8.12 Radian arc cotangent has two arccot func- 
tions. 

One is sign symmetric, but discontinuous at 0. 

The other is continuous, but not sign symmetric. 

That international standard is aimed at computer language standards bodies; so it may not meet your 
need for mathematicians. 

Fred J. Tydeman Tydeman Consulting 


The discussion continued in June but at the moment Ell refer to Johann's Titbits #32 pub- 
lished in this issue which deals with the problem - and we have too many other interesting 
requests which deserve to be presented. 
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Jan Vermevlen, Belgium 

Hello Josef, 

I have a question of one of my science colleagues. 

He has a set of some thousand quadruples (t,x,y,z) which show the position (x,y,z) of a parti- 
cle at a certain time t. 

Is it possible in Derive to produce a 3D graph of the path of this particle with an animation of 
the moving particle? Or should I use DPGraph? 

Thank you for giving me advice. 

Best regards, 

Jan 

Peter Schofield 

Hello Josef, 

It is nice to hear from you. Hope you are keeping well. I'm keeping myself busy teaching part-time at 
College. 

Re. your problem. I'm attaching a short Derive6 mth. file that might suggest a way of plotting a vector 
of 3D-points over a "time sequence". 

Open the File and set up a 3D-plot Slider Bar 1<= t <= 1000 with 999 intervals. 

Then aproximate and 3D-plot #3 or #4 using InsertPlot> Medium Points> Custom (Black Colours) 
(Might be better to turn Auto Change plot colours OFF as well). Takes about 60sec. on my laptop. 

The only drawback I can see is that the origin will also be plotted as well as the moving point, and so 
these plots look better with the Axes ON. I don't know how to get round this (minor) problem. 

Note that, in List3, 1 have stripped out the "time" coordinate of List4, since Derive can only plot 3D- 
points. However, there is still a notion of time steps in the successive points of List 3, and this is what 
#3 and #4 are using. 

All the best, 

Peter 

DNL 

Dear Peter, 

many thanks for your very valueable advice. It works. 

I attach a modified file - I don't need the 4-element vector and I defined two functions. 

The origin is nasty but it might happen that the range of the points does not contain the ori- 
gin, then it does not appear. 

Best regards 
Josef 
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File move j?ts . dfw 

This is the first attempt provided by Peter: 

( [x, X, X, x] 


#1: 1 ist4 := VECTOR 


#3: 


#4: 


-, x, -500, 500 


#2: list3 := list4- 


VECTOR C 
VECTOR C 


list3 


list3 


200 

0 0 0 
10 0 
0 10 
0 0 1 

■ IF(t = i, 1, 0), i, 1, 1000) 
■IFCt > i, 1, 0), i, 1, 1000) 


See the demonstration applied on 10 random points: 

I produce 10 randompoints in a 10 by 10 by 10 cube and 

trace the path using two functions derived from your algorithm: 

#5: VECTOR ([ [10 -RANDOM CD - 5, 10 -RANDOM CD - 5, 10- RANDOM (1) - 5]], i, 10) 


[-4.415, 

1 . 

,901, 

-1.113] 



[-4.415, 

1 . 

,901, 

-1.113] 

[-1.956, 

- 0 . 

,6297, 

-4.017] 



[-1.956, 

- 0 . 

,6297, 

-4.017] 

[-3.966, 

2. 

,134, 

-2.530] 



[-3.966, 

2. 

,134, 

-2.530] 

[-4.739, 

-3. 

,450, 

-2.002] 



[-4.739, 

-3. 

,450, 

-2.002] 

[-1.633, 

- 0 . 

,1096, 

1.452] 



[-1.633, 

- 0 . 

,1096, 

1.452] 





#7: 

pts := 





[3.636, 

-4. 

,949, 

1.404] 



[3.636, 

-4. 

,949, 

1.404] 

[3.425, 

-4. 

,347, 

3.263] 



[3.425, 

-4. 

,347, 

3.263] 

[0.3023, 

- 1 . 

,573, 

-0.5067] 



[0.3023, 

- 1 . 

,573, 

-0.5067] 

[-1.002, 


1.783, 

2.610] 



[-1.002, 


l . 733 , 

2.610] 

. [-3.930, 

-2. 

,955, 

-3.426] . 



. [-3.930, 

-2. 

,955, 

-3.426] . 


#3: 

movepts(list) := VECTOR(list -IF(t = i, 

i 

■ i — 

o 

■H 

DIM (] i st)) 

#9: 

tracepts(list) := VECT0R(1ist -IF(t > i 

i 

, 1, 0), i 

, DIM(]i'st)) 

#10: 

movepts(pts) 



#11: 

tracepts(pts) 
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David Siostrand 
Hi, 

I have spent some time thinking on Jan's problem having a point moving along a curve using 
a slider bar. In the attached file I have done it for a matrix containg 3 points. If tstep is con- 
stant you can do the same with matrices up to 1000 points but not more due to the limitation 
of the number of intervals in the slider bar. 

If you send me the actual matrix I can try to test my idea. 

Best regards, 

David 






" 1 -2 2 3 " 


Fi 1 e move_pts0 . dfw 


#1: 

A := 

2 3 4 5 


A is the matrix containing the 




.3 1 1 1 _ 


quadruples. If you plot #3 you 

re- 





ceive the path of the particle. 

in- 


fCt) 

:= VECTOR (A 

i 

sert a sliderbar with the below 


#2: 


t, i 


settings. Then plot #4 and then use 
the sliderbar to animate the move- 

#3: 

VECTORCfCt), t, 1, 

3) 

ment of the particle. 


#4: 

EHS 
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DNL: 

Your PPP-routine (Point - Path - Plotting) seems to work, that's great. I have produced a 
data set of a lot of random points and I 'll test your idea. 

Thank you on behalf of the DUG. (I thought to do something with a function - and you real- 
ized it in a very clear way!!, Congratulations). 

I was too enthusiastic about your idea and couldn't resist to make a try. I attach the file with 
a 100 point path using two functions based on your idea. 

It works!! 

File move_pts1 .dfw 


Dear David, 
Dear Jan, 



this is another version of „Voyage of the Point". 

The first column (“time column”) seems to be not necessary! 
Regards as ever, 

Josef 
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File move_j?ts2 . dfw 

#1: Notati'onDigits := 4 

mvClfst, t) := VECTORCIist , 1, 3) 

#2: t,i 

#3: VECTOR ([10 -RANDOM Cl) - 5, 10 -RANDOM Cl) - 5, 10 -RANDOM Cl) -5], i, 100) 

0.3547 4.269 -2.273 


-3.425 -3.422 3.532 


Result of #3 is assigned to pts2. 

Plot first the complete path by plotting pts2 point size small and points connected, then 
insert a slider bar with 1 <= t <= 100 (99 intervals) and plot with point size large 


In the meanwhile I received another mail from Belgium: 


Hello David and Josef, 

Sorry that I did not answer the last few weeks, but it has been (again) a bit too hectic in my neighbour- 
hood -©. 

Here is the data set that I received from my colleague. It is derived from a scientific experiment where 
a particle was sent through a whirling magnetic field (I think). 

Please note that, when you see the time sequence, there are big time intervals in between where the 
particle “disappears”. 

Can you test your solution with (parts) of the data? 


;j id a - -is* 


A2113 - f* 27234,1 



A 

B 

C 

D 


1 

time (ms) 

x (mm) 

y (mm) 

z (mm) 


2 

11773,4 

261,8 

34,2 

191,9 


3 

11782,3 

260 

44,9 

189,7 


4 

11787,6 

259,6 

54,2 

186,6 


5 

11791,8 

260,9 

61,1 

184,1 


6 

11797,4 

260,1 

72,4 

186,9 


7 

11802,8 

261,6 

81 

184,1 

8 

11806,8 

261,1 

87,5 

184,4 

9 

11811 

261 

94,6 

185,3 

10 

11814,6 

262,2 

101,1 

186,6 


7382 

188310,2 

270,9 

547 

218,6 

7383 

188315,4 

270,6 

550,2 

212,9 

7384 

188322,2 

271,2 

557,2 

216,8 

7385 

188330,3 

269,5 

563,4 

211,9 

7386 

188338,6 

268 

568,7 

209,8 

7387 

188349,5 

266,7 

572,2 

223,5 

7388 





7389 






h < ► n \Bladl/Blad2 / Blad3 / 


The problem kept Peter busy, too, improving his “Point Moving Procedure”: 
Hello Josef, 


I've been thinking some more about avoiding plotting the origin in these sliderbar applica- 
tions. 

Please find attached three Derive 6 files which go some way to solving this problem. 

You need "Approx before plotting" ON in both the 2D- and 3D-plot windows. 

Also, for the 2D-examples, use Points> Size (Large); Connect (Yes) 
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File movej?ts3 . dfw 


In 2D: 100 randompoints in an S by S square and trace the path using four functions 
derived from the improved algorithm: 

tnoveptsOist) := VECTOR (list -IF(t = i, 1 , co), i, DIM(list)) 


#1: i 


tracepts(list) := VECTOR(list -IF(t > i, 1, co), i, DIM(list)) 
#2: i 


#3: pts :i VECTORCLLS-RANDOMCl) - 4, S-RANDOMCL) - 4]], i, 100) 


#4: 


ptsl :i VECTOR C 


[0, 0], pts 


i , 1 J - 


, i, 1, DIMCpts)) 


Set up a 2D-sliderbar: l<=t<=100 with 99 intervals, 

I plot the point list (grey and size Medium), then: 



movepts(pts) 

Shows one point after the other by moving 
the slider. 


-5 


tracepts(pts) 

Shows the first t points by moving the slider. 




movepts(ptsl) and tracepts(pts2) show the points + pointers to them. 
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File move_pts4 . dfw 


#1: 

movepts(list) := VECTOR(list -IF(t = i, 1, «0 , i, 

i 

DIM (list)) 

#2: 

tracepts(list) := VECTOR(list -IFCt > i, 1, m) , i 

i 

, DIM (list)) 

#3: 

movepts(pts) 


#4: 

tracepts(pts) 



In 2D: SO points on M y = x A 3/10 M containing a data point at the origin. 
How to adapt the data list so that this point is not plotted 
each time when using a slider bar? 


|nudge(v) := 

If NUMBER_TYPE?(v) 

If v = 0 
10 A C-10) 
v 

VECTQR(nudgeCvii) , 


DIM(v))| 


# 6 : 


pts 


f 


\~ VECTOR 

V 


3 

i - 40 Ci - 40) 
10 ’ 10000 




SO 


Plot #6 as grey mid sized not conntected points. 
Set up a 2C^sl iderbar : l<=t<=S0 with 79 intervals 
and observe the point "Walking the Line 11 . 


#7: movepts(nudge(pts)) 



#3: tracepts(nudge(pts)) 



Comment: Tracing the points only on a pa- 
rameter defined curve is easier done by: 

and applying a sliderbar for i. 

I wanted to make the “path” of the point visi- 
ble and experimented with Peter's functions. 


3 

i - 40 Ci - 40 ) 
10 ’ 10000 


See my results in 

tracepts_conn (nudge (pts)) 


Josef 
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#9: 


tracepts_conn(list) := VECTOR 


(list -IF(t > 

i 


1, m)) , (list 


■ IF(t >i+l, 1, m3) 


i + 1 




DIM(list) - i) 


#10: tracepts_conn(nudge(pts)} 


File move_pts5 . dfw 



100 randompoints in a 10 by 10 by 10 cube and trace the path using four functions 
derived from the improved algorithm: 

movepts(list) := VECTOR(list -IF(t = i, 1, cc), i, DIM (list)) 

#1: i 

tracepts(list) := VECTORClist -IF(t > i, 1, co) , i, DIM C list)) 


#2: i 

#3: pts := VECTOR ([[10 -RANDOM (1) - 5, 10- RANDOM (1) - 5, 10- RANDOM (1) -5]], i, 100) 


ptsl := VECTOR ( 


#4: 


[0, 0, 0], pts 




1,1, DIM (pts)) 


Set up a BC^sl iderbar : l<=t<=100 with 99 intervals. 

#5 : pts 

#6: 


novepts(pts) 


#3: 


tracepts(pts) 







D 
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#10: [tracepts(pts) , tracepts(ptsl) ] 



Now it is time to return to the given problem: tracing the particle through a magnetic field. 
First of all I transferred the first 1000 positions from Excel via Word to DERIVE, taking col- 
umns 2, 3 and 4. (See DNL# how this can be performed. Another contribution (next DNL) 
will explain how to assign a data matrix to a variable using a DERIVE program.) 

File gegevens . dfw 

#1: NotationDigits := 4 

mv(list, t) := VECTOR Cl ist , i p 3) 

#2 : t , i 

#3: gegevl := [[261.8, 34.2, 191.9], [260, 44.9, 189.7], [259.6, 54.2, 186.6], [260.9, 61.1, 184.1], 

[260.1, 72.4, 186.9], [261.6, 81, 184.1], [261.1, 87.5, 184.4], [261, 94.6, 185.3], [262.2, 101.1, 
186.6], [262.1, 107.1, 179.5], [261.8, 113.5, 180.4], [263.6, 118.6, 182.3], [262, 121.5, 183.2], 


This is to find appropriate values for the bounding box : 

#4: VECTOR ([MIN (gegevl 1 1 k) , MAX(gegevluk) ] ■ k, 3) 


# 5 : 


241.3 266 

34.2 572.3 

_ 174 216.7 _ 
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next produces the path of the particle (1st 1000 positions) 

# 6 : 


gegevl 



#7 : BBBBBiHSH 

Introduce a slider bar for t and the plot #7, 

Highlight #7 and Insert > Plot > point size Medium or Large, 

Colors Custom (I set all colors RED). 

Voila, it works! ! 

We see the particle on its position #72 and then 

just before finishing the first part of its journey through the magnetic field. 





D 
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#9: 

#10: 

# 11 : 


pa.thstep3d(list, t) := 
mv_e("h'st, t) := VECTOR (list 


VECTOR (list , i, 3), VECTOR ("list , i , 3)' 

t , i t + 1 , i 


i. 3) 


t + l,i 

[pathstep3d(gegevl) , mv_e(gegevl)] 



A Random Walk 

I was so excited about our success that I wanted to produce a simulation of a simple random 
walk in the 2D-plane: 


#1: 

mv2dClist 1 t) := VECTORflist , i, 2) 

t, 1 





#2: 

mv3d(list, t) := VECTOR(list , i s 3) 

t, i 





#3: 

pathstep2d(list, t) := rVECTOR(list , i , 

L t, i 

2), 

VECTOR ("list 

t 

+ 1,1 

i, 2) 

#4: 

pathstep3dClist, t) := rVECTORflist , i, 

L t 3 i 

3), 

VECTORCIist 

t 

+ I,/ 

UJ 

W 


randwalk(n, m, dummy, p, i, k) := 






Prog 

m := [ 1 , 0; -1, 0; 0, 1; 0, -1] 

dummy := RANDOM CO) 

P := [0, 0] 
trace := [p] 
i := 1 

#5 : Loop 

If i > n 

RETURN Fl result in matrix trace' 1 
k := RANDOM (4) + 1 
p := p + mj,k 

trace := APPENDCtrace, [p]) 
i : + 1 

#6: randwalk(1000) = result in matrix trace 

#7: trace 


#3: 


[nv2d(trace)| 
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#6: randwalk(lGOO) = result in matrix trace 

#7: trace 

#8: mv2d(trace) 



#9: randwalk(lOOO) = result in matrix trace 

#10: trace 

#11 : p at h s t e p 2 d ( t r ac e ) 



We can perform the walk on the Tl, too. Use the following program for observing a point ran- 
domly strolling around on the TI-92 or Voyage graph screen: 



D 
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randomwal k(n) 

Prgm 

Local di rs , i , j , p , pn 

FnOff :ClrDraw 

[[1 ,0] [-1 ,0] [0,1 ] [0, "I ] ]-di rs 

“23 . 8^xm1 n : 23 . 8->xmax 

“10. 2->ym1 n : 1 0 . 2->ymax 

[ [0,0] ]-p: [ [0,0] ]-pn 

Circle p[1 ,1] ,p[1 ,2] ,0.5 

©Pause 

For 1 , 1 , n 

Circle p[1 ,1] ,p[1 ,2] ,0.5,0 

©Line p[1 ,1] ,p[1 ,2] , pn [ 1 ,1] , pn [ 1 ,2] 

p+dl rs [rand (4) ]^pn 

Line p[1 ,1] ,p[1 ,2] , pn [ 1 ,1] , pn [ 1 ,2] 

For j , 1 , 5 

Circle pn [1 , 1 ] , pn [1 , 2] , 0. 5 

EndFor 

©Pause 

pn+p 

EndFor 


r'FiTTiriY fet y fjt y fht y fe y FfiT T "i 
|fllgebra|Calc|0ther|PrgmI0|ciean Up| 1 

■RandSeed 1234 

□one 

■ randual k(2G0) 

Done 

r an dual k < 2 □□ >1 

MAIN RAD EXACT 

FUNC E/30 




A fine application of a Random Walk is the Brownian Motion. A Random Walk-hypothesis models the 
share prices based on a Brownian Motion with expected value and volatility of the yield of the share as 
parameters. 

See three simulations of 800 steps how a share with a price of 50 will develop (mean and volatility of 
yield = 0.01 and 0.4) 

(Part of a possible future contribution, Josef) 


geombrown(50, 800, 0.01, 0.4) 
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I received an interesting mail from Heinz Rainer Geyer concerning a random experiment 
which I had used very often as an introductory example for probability theory: I brought three 
dice with me in class and asked the students what they would like to pay for one game: roll- 
ing the dice and win 1 through 6 EURO for a sum from 13 to 18 respectively. 

They all made their offers . . .. 

Then we made a number of experiments, noted the incomes and the payments of the bank- 
holder and finally we calculated the expected value of the payment. What would you like to 
bet? For calculating the mean we needed the answer to Heinz Rainer's question: 

Heinz Rainer Geyer, Germany 

Hallo lieber Josef, 

vor einer Stunde wollte ich dich noch mit einem Problem belastigen, das ich jetzt zumindest ansatz- 
weise gelost habe. Es geht um die alte Aufgabe der Verteilung der Augensummen bei n gleichzeitigen 
Wiirfeln. 

Ich habe lange keine eindimensionale Liste aller Ergebnisse bei 3 Wiirfeln aufbauen konnen. 

Die jetzige Losung ist auch nicht gerade elegant und schreit formlich nach einem anderen Aufbau. 

Hast du dafur einen Tipp? 

Das zweite Problem ware eine Verallgemeinerung der Strukturen, so dass man nur noch die Anzahl 
der Wiirfeln als Variable verwenden muss, um die Verteilung der Augensummen zu erhalten. 

Kennst du eigentlich einen Algorithmus zur direkten Berechnung der Haufigkeiten der Augensummen, 
zumindest bei 3 Wiirfeln? Ich habe bisher nichts gefunden. 

(Heinz Rainer wants to automate the well known problem to find the frequency distribution of 
the sum of three (or more) rolled dice. He had some problems to build a list of all possible 
cases for three dice and asks for generalisation for any number of dice rolled. 

Finally he asks if I knew any algorithm to obtain directly the distribution of the sums. See the 
last lines of his DERIVE file. 

#19: wurf4 := APPEND(VECTOR(APPEND(VECTOR(APPEND(VECTOR (VECTOR ([t, u , v, w], w, 1, 6), v. 1, 6)), u, 1, 6)). 

t, 1, 6)) 


(DIMCvO 


#20: 

zaehl_4(n) 

:: DIM 

SELECT 


I w_ = 

n, w_, wurf4 
















i=l i 

JJ 










#21: VECTOR ([n, 

zaeh"l_4(r0], 

n, 

4, 24)' 












'45 6 

7 

8 

9 

10 

11 12 

13 14 15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

#22: 

. 1 4 10 

20 

35 

56 

80 

104 125 

140 146 140 

125 

104 

80 

56 

35 

20 

10 

4 

1 


DNL: 

Lieber Rainer, 

ich schicke Dir ein kleines Programm, das fur eine allgemeine Zahl von Wiirfeln alle mogli- 
chen Wiirfe mit den Augensummen und deren Haufigkeiten angibt. 

Ich bin (fast) sicher, dass es noch eleganter geht. 


D 
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I sent a short program which returns for any number of dice the requested distribution, but I 
was quite sure, that it must be able to find a more elegant way to tackle the problem. 

Heinz Rainer: 

Hallo lieber Josef, 

das ist ein interessantes Programm, da muss ich mich mal einarbeiten. Mit C0DES_T0_NAME() usw. 
hab ich mich bisher noch nicht befasst. 

Offensichtlich hast du mit APPEND(1 i st , [wurf ]) deine Struktur zusammensetzen konnen. 

Ich wollte ja aus der Ergebnismenge mit 6 n Elementen die Summen mit SELECT abzahlen. 

Dazu miisste ich zu jedem Element des vorherigen Vectors [wl,w2,...,w n _i] die 6 neuen Wiirfe anfu- 
gen. 

Irgendwie vermisst man dabei die FOR Schleifen, aber mit LOOP sollte das ja auch gehen. 

Heinz Rainer worked through the program and noted that he had not worked with 
C0DES_T0_NAME() before. He missed the FOR-loops but LOOP should work in a similar way. 

Heinz Rainer: 

Heureka, 

Mann, das lies mich jetzt nicht mehr los! Es hat mich ein ganzes Wochenende und den Nachmittag 
gekostet, aber es ist der direkte Zugang, den ich mir vorgestellt hatte. 

Elegant ist das sicher nicht. 

Hast du das Problem schon an Johann gemailt? Er macht bestimmt einen rekursiven 4-Zeiler draus? 

Liebe ( ini Be 
Rainer 

Danke fur deine neue Datei, ich habe sie in der letzten Zeile gerade ein wenig umgestellt, um die Hau- 
figkeiten auch zu plotten. So sieht das naturlich viel kompakter aus als meine vielen LOOPs. 

Finally Heinz Rainer worked hard a full weekend to find a direct approach but he admitted, 
that his solution would not be very elegant. Then he asked: “Did you send the problem to 
Johann (Wiesenbauer = Mr. Titbits). I am sure that he will answer with a recursive 4 line 
function”. 

Great idea, I asked Johann and received an answer within a few hours. 

Lieber Rainer, 

und hier ist er, der "Einzeiler" von Johann Wiesenbauer. 

Er hat mir diese mail geschrieben, die Funktion dazu zu basteln, war nicht schwierig. 

Liebe GriiHe und viel Spall beim Zocken, 

Josef 

This is Johann's answer: it is not 4 lines but only one DERIVE function used in a clever way 
to give the distribution of the sum of three dice. It was not difficult to generalise this function 
for any number n of dice. 
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Von: Johann Wiesenbauer [ mail to : j . wiesenbauer6tuwien . ac . at ] 

Lieber Josef, 

Was Deine Anfrage, so ist die tatsachlich ganz einfach zu beantwor- 
ten . 

Simplif iziere mal 

table (poly_coef f ( (x+x A 2+x A 3+x A 4+x A 5+x A 6) A 3 ,x ,k) ,k,3,18) 

was eine Tabelle fur die absoluten Haufigkeiten der Augensummen bei 
3 Wurfeln liefern sollte. Du kannst Dir dann leicht selbst uberle- 
gen, warum das funktioniert und wie man das allgemein auf n Wurfeln 
verallgemeinern kann . 

Herzliche GriiDe, 

Johann 


2 3 4 5 6 4 

#1: TABLE CPOLY_COEFFCCx + x + x + x + x + x ) , x, k), k, 4, 24)' 


#2: 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 


. 1 

4 

10 

20 

35 

56 

80 

104 

125 

140 

146 

140 

125 

104 

80 

56 

35 

20 

10 

4 

1 


And here is the function of your desire (works in an instant): 

2 3 4 5 6 n 

#3: hwsumme(n) := TABLE (POLY_COEFF((x + x +x +x +x + x ) , x, k), k, n, 6-n)’ 



4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

hwsumme(4) = 

. 1 

4 

10 

20 

35 

56 

80 

104 

125 

140 

146 

140 

125 

104 

80 

56 

35 

20 

10 

4 

1 


You see Johann's function as core of a TABLE. It is a wonderful tas for students to investi- 
gate this function and find out why it works!! 


Finally I'd like to present my - DERIVE program - based on counting the frequencies in its 
last version. It is nice to plot the transposed results with points connected 
(file dice_sum . dfw). 

Additionally I wanted to solve the problem for the Tl, too. So I invented the poly_coeff- 
function for the Tl. I find that this might be an interesting problem for students. Applying this 
function it is not difficult to find a function for returning the distribution of the sums 

dice_cnt (n) . 

(For those of you who don't now poly_coeff: poly_coeff (u,v,k) returns the coeffi- 
cient of v k in the polynomial u.) 

Lieber Rainer, 

ich bin eben dabei, den nachsten DNL zu "komponieren". Bei dieser Gelegenheit habe ich 
mein Wurfelprogramm noch etwas abkurzen konnen, da ich das Zahlen der Haufigkeiten 
wesentlich einfacher gestaltet habe, ohne die Grundidee des Programms zu andern. 

Es ist nett, die transponierten Ergebnismatrizen "verbunden" zu plotten. 

Johanns Trick mit poly_coeff habe ich fur den V200/TI-92 programmiert. 
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di'ce_sh(n, start, end, wurf, k, wcode, hlist) := 

Prog 

hlist := VECTORCO, ], 1, 6-n) 

start := C0DES_TOJ^AME (VECTOR (49, i, n)) 

end := CQDES_TO_NAME (VECTOR (54 , i, n)) 

wurf := start 

Loop 

If wurf > end exit 

wcode := NAME_TO_CODES(wurf) 

#22: k := 1 

Loop 

If wcode | k < 49 v wcode. |k > 54 exit 
k :+ 1 
If k > n 

h 1 1 s t i (K NAM E_TG_CG DE S (wu r f ) ) - n-43) : + 1 
If k > n exit 
wurf : + 1 
k := n 

RETURN [VECTOR (i , i, n, 6-n), hlist|[n 6-n]] 


#23: 

dice_sh(4) 





4 5 

6 

7 

8 

9 

#24: 

. 1 4 

10 

20 

35 

56 

#25: 

dice_sh(5) 





5 6 

7 

8 

9 

10 

#26: 

. 1 5 

15 

35 

70 

126 


27 

28 

29 

30 



35 

15 

5 

1 . 



The distribution of the absolute frequencies of the sum of three, four and five dice 



^ YFStYFHtV fe YFfiTT 

^ | C ■=■ n t r ■=■ 1 1 1 Q | J a r| F i n d . . . | N ■=■ d e | 


: poly_cf Cu ? x, rO 
: x ? n > I x=G 


t j— I C ■=■ n t r o 1 1 1 0 1 U a r | F i n d . . . | N o d e | 


: dice_cnt(rO 
: Func 
: Local pis ? f req 
: seq(! k ? k ? n ? 6*n>-*pt-£ 

: seq(poly_cf C Cx+x^2+x A 3+x A 4+x A 5+x A 6^n ? x 


_ . pol' 

, k> , k r n, 6*n>+freq 
: 1 ist> mat (augment- Cpts ? f req) , 5*n+l ) 
: EndFunc 




24 


■ poly_cf((b 3 + 2 ■ b 2 + 2) 3 , b ? 4) 

■ expand((b 3 + 2 b 2 + 2) 3 ) 

b 9 + 6 b 8 + 12 b 7 + 14 b 6 + 24 b 5 + 24 b 4 ► 


1 poly. 




a-x° + b*x 


:) ? x ? 7] 


5 ■ a ■ b 


expand< <a*x A 3+b*x> A 5 > 


FITTHT fT? Y FSt Y FHt y f? y f_ 
f — | R 1 g e b r a | C a 1 c 1 0 1 h e r \ F r g n 1 0 | C 1 e a n IJ p 


f'[ ( a ■ x 3 + b ■ x ] ^ , x , 7 ] 


D 


5 ■ a ■ b 


1 poly_cfl 
1 dice_cnt(2) 

f2 3 4 5 6 7 8 9 IQ 11 12' 


E 


23456543 
1 dice_cnt(4) 
f4 5 6 7 


8 


1G 11 


2 

12 13 


n 


Ll 4 10 20 35 56 8Q 1Q4 125 14Q 

MlilN FiliD EKflCT FUMC E/30 
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Titbits from Algebra and Number Theory (32) 

by Johann Wiesenbauer, Vienna 

As you may know (or read in this issue), there has been a quite interesting discussion on 
the Derive newsgroup recently. It was started by Jim FitzSimons, who wanted to know 
how Derive does find a certain algebraic expression for cos(7i/17). Well, as I pointed 
out in my answer, I had already dealt with a related question in my "Titbits" (in the DNL 
#20), and later I even repeated the main points of that derivation in an attached De- 
rive-file, but the further discussion made it desirable to look into this topic again from 
a more general (you could also say more algebraic !) point of view. In particular, when 
Albert Rich asked for a radical expression for cos(7t/257), it was clear for me that I 
should rewrite that article and have a try at this new challenge. What you read in the 
following is the outcome of this plan, which is neither a sweeping success nor a total 
failure though as to this goal. 

In order to make everything as transparent as possible, I'll give all basic formulas for a 

general Fermat prime p, i.e. a prime of the form 2 2 " +1 for some natural number n. On 
the other hand, as for examples, to save space I'll consider here only the case p = 17 
rather than p = 257, which is dealt with though in the Derive-file accompanying this 
column. 

To begin with, we need a so-called primitive root g mod p, i.e. an elementg, whose 
residue class gmodp is a generator of the prime residue class group Z p *. (Note that 
such an element always exists according to a theorem by Gauss!) In Derive the smallest 
g with this property can be obtained like this: 


PRIMITIVE_R00T (p) 


g := 3 


Next we form a 1-class partition aO containing all odd numbers below p. 


aO := VECTOR IF (ODD? (O, a^, p - O, a^, ITERATES MODCg-x , pj , x, 1, 


aO := [[1, 3, 9, 7, 13, 5, 15, 11]] 


What a strange way to compute those odd numbers, you might be saying! Of course, you 
are right, but the secret is the very specific order of these elements as will become 
clear later. 
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Furthermore, we will need the vector bO of all sums of the terms 2cos(k7i/p) for all k 
contained in some list of a . Note that when referring to 2cos(k7i/p) we will use the 
term z k -z p_k in the following, where z = e 1,t/p or z is rather a root of the irreducible 
polynomial (z p +l)/(z + l)over Q. (Numerical comparisons below will make sure that the 
"right" root is taken from the p different roots!) 

k p - k 

bo := VECTOR (I(z - z , k, O, a_, aO) 


16 15 14 13 12 

z +z -z +z -z 


11 10 9 

z - z + z 


z + z - z + z 


Furthermore we will need the vector c of the corresponding real numbers represented 
by those sums. (Note that we deal in the following only with approximative values of c , 
which is sufficient for our purposes!) 



f 

' k-7T ^ 

> 

> 

2-1 

COS 

^ P J 

, k, a_ 
> 

, a_, aO 

> 


APPROX(cO) = [1] 

Finally, we introduce a vector dO that contains the "exact" values of the components of 
cO . At the beginning, dO is always [1] as suggested by the approximative value of cO 
above. 

dO := [1] 

Before going on, we write a routine that does all the steps above automatically, if the 
parameter p is initialized by p. 

init(p) := 

Prog 

p := p 

g := PRIMITIVE_ROOT (p) 

aO := [VECTOR (IF (ODD? CO, a_, p - O, a_, ITERATES (M0D(g -x , p) , x, 1, (p - 3)/2))] 
bO VECTOR(I(z A k - z A (p - k) , k, a_) , a_, aO) 
cO := VECTOR (2. X(C0S(k.7T/p), k, O, a_, aO) 
dO := [1] 


Note that all variables in init(p) (except for p) are global variables, which can be viewed 
and also used after the call of the routine. (Actually one of the rare cases, where the 
use of global variables makes sense in Derive!) 

We now divide each list in aO into two equally sized parts by taking the oddindexed 
numbers and the evenindexed numbers thus getting a new list a . 


a := APPEND(VECTOR([" a_ 


DIHCO] [2, 4, 


DIM (a 


tO]T 


a_, aOJJ 


11 
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Updating bO and cO as well yields the following new values. 



f 

r 

k 

P - k 


p ) 

z + 1 


b := VECTOR 

REMAINDER 

ICz 

V 

- Z , 1* 

i, o, 

z + 1 J 

, a_, a 

j 


I" 14 12 11 10 7 6 5 3 14 12 11 10 7 6 5 3 

b Lz +z -z + z -z+z-z-z+l.-z -z +z -z +z-z+z+z 


VECTOR 2-1 COS 


APPROX(c) = [-1.561552S12, 2.561552S12] 

Now consider two adjacent elements u and v of b such that u = v 2i _! and v = b 2i for 
some i = l,2,...,dim(b)/2. In order to compute d, i.e. the vector with the exact values, 
we need the crucial fact that both u + v and u*v, when reduced mod (z p +l)/(z + l) are 
linear combinations over Q of certain elements of bO . (Actually even with integer co- 
efficients, but we don't make use of this!) 

This is trivial as far as u + v is concerned. (Remember, we got the elements of b by 
"splitting up" the elements of bO into "two halves".) As for the product u*v, let’s 
check it here for the only pair of b . 


= [-4] 


In other words, the numerical values of u and v can be obtained as solution of some 
quadratic equation, e.g. in the case at issue: 



r p 

> 


z + 1 

DIM(b) 

REMAINDER 

b -b , 

. i. 1. 


^ 2 - i - 1 2-i z + l> 

2 J 


FIRST (SOLUTIONS ([u + v = 1, u-v = -4], [u, v])) 


yl7 1 1 Vl? 

2 2 ' 2 2 


Actually there are always two solutions to a system of equations u + v = r, uv = s with 
known values of r and s , which only differ by the order though, and to find the "cor- 
rect order" we must resort to the numerical values stored in c . 


solution( 
If cl > 

[<V( 


r, s, cl, c2) 


4-s) + r)/2, (r - ^(r A 2 - 4-s))/2] 


[(r - V(r A 2 - 4 ■ s))/2 , (V(r A 2 - 4- 
d := APPENDfvECTORfsolutiona, -4, c 


r)/2] 


DIH(c) 


), i, 1, 


1 Vl7 Vl7 


2-1 - 1 2-1 
Jll 1 


Before going on, let's save the values of a,b,c,d to a0,b0,c0,d0, respectively, and com 
bine all the assignments above in single routine called nextstep( ): 
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[aO a, bO b, cO c, dO d] 

nextstep(show := true, s_) := 

Prog 

a := APPEND(VECTOR([a_4,[l, 3 DIM CO], a_4,[2, 4 DIMCOJ], a_, aO)) 

b := VECTOR (I(z A k - z A (p - k) , k, O, a_, a) 
c := VECTOR (2- I(C0S(k.7T/p), k, O, a_, a) 

s := VECTOR (REMAINDER (bj. (2-i - 1) .bj,(2. i) , (z A p + l)/(z + 1)), i, 1, DIM(bO)) 

s_ (z A p + l)/(z + 1) ■ SUBST (s_, z, 0) 

s_ := VECTOR (VECTOR (P0LY_C0EFF(t_, z, O, a_, aO COL 1) , t_, s_) 
d APPEND(VECTOR(solution(dOik, s_4,k.d0, cj,(2.k-l), cj,(2.k)), k, 1, DIM(bO))) 

[aO := a, bO := b, cO := c, dO := d] 

If show exit 
FIRST (d) /2 

The only "tricky" part is to represent the elements of s_ (in the program after the very 
first assignment) as linear combinations of (z p +l)/(z + l) and the "old" basis bO. This is 
done in the following two assignments, after which s_ contains the vectors with the in- 
teger coefficients of all those linear combination with respect to bO. Exactly these 
coeff icients are needed to obtain the new "exact" vector d from dO . As an output af- 
ter each step (if wanted, i.e. show is true), we get dj/2, which is the sum of all 
eypressions cos(k7i/p), where k ranges in that class of the partition of a , which con- 
tains 1. At the end, this class contains the number 1 only and hence we get the result, 
we are longing for, namely cos(7i/p) . 

Let's check it by continuing using this routine, which takes only tenths of a second on a 
fast PC. 

nextstepO 

7(34 - 2-717) 7l7 1 

S S S 

nextstepO 

7(7(3S-7l7 + 170) + 3.717 + 17) 7(34 - 2-717) 7l7 1 

+ - + 

S 16 16 16 

In fact, Derive acknowledges the following equality (in approximate mode though!) 


IT 1 ) 7(7(38 .717 + 170) + 3.717 + 17) 
■ - 

7(34 - 2-717) 

717 


k L 17 J S 

16 

16 

kD 

1 

h 


Now, what about the case p=257? Again the computations are very fast (though taking 
seconds rather than splits of a second), but as Jim found out independently, Derive 
can't cope with the huge expressions short before arriving at the final expression. 
Frankly, if you look at the huge expressions in the Derive-file you really wonder why 
Derive didn't throw in the towel earlier! 

It must be said, that Jim FitzSimmons obtained a result, too. He wrote in his last email: “Johann, I can 
not get Derive to do the last step. The expression is too long for Derive to handle. The example is 
attached. (cos257b.mth) 

Regards, Jim“ 
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See the last lines of Mr Titbits' file deriving COS(ti/257): 


#29: nextstepO 

2 

VCdl + 4-Cdl7 + d30)) + dl 

#30: 

4 

The next two assignments (the sense of which will become clear only later!) must be again simplified before 
going on, although I have delete the putput to save space. 


ddl 

:= SUBSTCd , 

, ds, 

d_) 

#31: 

1 


dd8 

:= SUBSTCd , 

, ds, 

d_) 

#32: 

8 



For the last step we must repeat the whole procedure with different notations. 

#33: dds := VECTOR (APPEND(dd, STRING (k)) , k, 1, DIM(d)) 

#34: dds := [ddl, dd2, dd3, dd4, dd5, dd6, dd7, ddS, dd9, ddlO, ddll, ddl2, ddl3, ddl4, ddl5, ddl6, ddl7, 

ddlS, ddl9, dd20, dd21, dd22, dd23, dd24, dd25, dd26, dd27, dd28, dd29, dd30, dd31, dd32, dd33, 

dd34, dd35 , dd36, dd37, dd38, dd39, dd40, dd41, dd42, dd43, dd44, dd45, dd46, dd47, dd48, dd49, 

dd50, dd51, dd52, dd53, dd54, dd55, dd56, dd57, dd58, dd59, dd60, dd61, dd62, dd63, dd64] 


#35: dO :i dds 

#36: nextstepO 


2 


VCddl + 4 -dd8) + ddl 


Hence this is the final symbolical formula for cos(pi/257). (Note the difference between "ddr',"dd8" and ddl, 
dd8, respectively, which is not shown on the screen!) 

r 2 


#38: cos_pi_over_257 :z SUBST 


VCddl + 4 -dd8) + ddl 


[ddl, dd8] , [ddl, dd8] 


Unfortunately Derive can't carry out this substitution due the size of the expressions involved. (If you are a 
perfectionist though, you can still go one step further and substitute "dd8" by dd8 alone and leave "ddl" 
unsubstituted, but I left it at that. ) Well, it's like a mountain climber must feel, if he is 100 m below the peak 
of the Mount Everest, but can't reach it for some reason, isn't it? 


At least, we can see that the result is apparently correct. In fact, when approximating this term with an 
increased precision of 100 digits we get 

#39: APPROXCSTRING (cos_p i _ove r_25 7) , 100) 

#40: 0.9999252866697325552128053591455562406916298296429212011881377099710071475277558301156004129496104709 


On the other hand, we have 


r r 

r n Y| 

"N 

I STRING COS 

l 257 JJ 

, 100 

J 


#42: 0.9999252866697325552128053591455562406916298296429212011881377099710071475277558301156004129496103685 


which shows a remarkable coincidence and is sort of a "proof" of our assertion. 


